Phase separation in star polymer-colloid mixtures 
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We examine the demixing transition in star polymer-colloid mixtures for star arm numbers 
/ = 2, 6, 16, 32 and different star-colloid size ratios. Theoretically, we solve the thermodynami- 
cally self-consistent Rogers- Young integral equations for binary mixtures using three effective pair 
potentials obtained from direct molecular computer simulations. The numerical results show a spin- 
odal instability. The demixing binodals are approximately calculated, and found to be consistent 
with experimental observations. 
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The study of mixtures of hard colloidal particles and 
nonadsorbing polymer chains has received a great deal of 
recent attention, both experimentally and theoreti- 
cally |}]-|3]]. However, very little is known about the phase 
behavior of the much richer mixture of colloidal particles 
and star polymers, the latter being macromolecules con- 
sisting of / polymeric chains covalently attached on a 
common center. The theoretical approaches to the study 
of colloid-polymer mixtures were largely based on the 
Asakura-Oosawa (AO) model, in which the chains are 
envisaged as noninteracting spheres experiencing a hard- 
sphere repulsion with the colloids. Subsequently, the 
system can be mapped onto an effective one-component 
fluid featuring the so-called depletion interaction between 
the colloids, mediated by the ideal chains Though 
the AO model provides an excellent benchmark for such 
systems, recent theoretical studies [^,^| and comparisons 
with experiments Q indicate that the assumption of non- 
interacting chains leads to quantitative discrepancies be- 
tween the two. Hence, a systematic effort to derive more 
realistic chain-chain j?| as well as chain-colloid |||| in- 
teractions has already been undertaken. 

For the case of a star polymer-colloid mixture, the 
need to employ realistic interactions is even more ap- 
parent. Star polymers with a large functionality / are 
stiff particles and their mutual interactions cannot 

be ignored, even as a first approximation. Indeed, / is 
the parameter which governs the softness of the star-star 
as well as of the star-colloid interaction, as we demon- 
strate below, and which provides the natural bridge be- 
tween the colloid-chain mixtures (at / = 1 or 2) and the 
binary hard sphere mixture, formally equivalent to the 
limit f — ► oo. Moreover, by changing / or the degree 
of polymerization N of the star arms, a large domain 
of size ratios between the components can be covered. 
When the stars are of comparable size with the colloids 
a depletion picture of the mixture will be less accurate 



than a full, two-component description because the for- 
mer ignores the effects of many-body terms. In this pa- 
per, we propose analytical expressions for the star-star 
and star-colloid interactions for all / values and for a 
large range of size ratios between the two. These are 
based on the one hand on accurate, monomer-resolved 
computer simulations and on the other on theoretical ar- 
guments regarding the functional form of these interac- 
tions. Using these expressions, we make theoretical pre- 
dictions on the mixing-demixing (or 'gas-liquid') tran- 
sition in star polymer-colloid mixtures. Our measured 
demixing curves are in good agreement with the theoret- 
ical predictions. 

We consider a binary system with N c colloidal spheres 
of diameter a c (radius R c ) and N s star polymers, charac- 
terized by a diameter of gyration er g (radius of gyration 
R g ) and an arm number /. For the special case / = 2, 
star polymers reduce to linear polymers and we reach 
the limit of colloid-polymer mixtures [|lM@|. The total 
particle number is N = N c + N s . Let q — cr g /a c be the 
size ratio and rj c — -y-^ cr 3 and r/ s = a 3 the packing 
fractions of the colloids and stars respectively, enclosed 
in a volume V. 

Experimentally, we studied two sets of star polymer- 
colloid mixtures consisting of poly(methylmcthacrylate) 
(PMMA) particles and poly(butadiene) (PB) star poly- 
mers with size ratios q ~ 0.49 and q « 0.18, respec- 
tively. PMMA particles were synthesized following a 
standard procedure O). Stock suspensions were pre- 
pared either in cis-decahydronaphthalene (cis-decalin) or 
cis-decalin/tetrahydronaphthalene (tetralin) mixture as 
an index-matched solvent. These systems have been es- 
tablished as hard sphere models [Q. The volume frac- 
tion rjc was calibrated using the onset of the hard sphere 
freezing transition, taken to be at r\ c — 0.494 and ob- 
served as the nucleation of iridescent colloidal crystals. 
The PB star polymers were prepared by anionic poly- 
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merization following an established procedure |12|,|13[. 
Star arms were synthesized by polymerizing butadiene 
with secondary butyl lithium as initiator. The resulting 
living polymer chains were coupled to the chlorosilane 
linking agent having ideally 6, 16 and 32 Si-Cl-groups. 
The molecular weights M w of the PB arms were ad- 
justed to give star polymers with values of (i? 2 .) 1 / 2 = 
0.0172M°- 609 /" - 403 @ as close to 50 nm as possible. 
A linear PB polymer (/ = 2) was prepared as a refer- 
ence system. The particles and star polymers were char- 
acterized using light and small angle neutron scattering 
(SANS) Q. The results are summarized in Table |[ 

Samples were prepared by mixing PMMA suspensions 
with PB stock solutions. Each sample was homogenized 
by prolonged tumbling and allowed to equilibrate and ob- 
served by eye at room temperature T — 25 °C In 
all samples with q w 0.49, addition of polymer to sus- 
pensions with r] c ~ 0.1 — 0.4 brought about, successively, 
phase separation into colloidal gas and liquid (or demix- 
ing), triple coexistence of gas, liquid and crystal, and 
gas-crystal coexistence. In samples with q 0.18, ad- 
dition of polymer first led to fluid-crystal coexistence; a 
metastable gas-liquid binodal buried inside the equilib- 
rium fluid-crystal coexistence region was encountered at 
higher polymer concentrations [T^ ]. In all cases, demix- 
ing started within several hours, crystallization within 
two days. Here we focus on the demixing transition. 

Theoretically we model the thermodynamics of the 
mixtures on the level of pair potentials between the 
two mesoscopic components, having integrated out the 
monomer and solvent degrees of freedom. Thus, three 
pair potentials are used as inputs for thermodynamically 
self-consistent integral equations which are closed with 
the Rogers- Young (RY) scheme p8| . The colloid- colloid 
interaction at center-to-center distance r is taken to be 
that of hard spheres (HS): 

The effective interaction between two stars in a good 
solvent was recently derived by theoretical scaling ar- 
guments and verified by neutron scattering and molec- 
ular simulation, where the monomers were explicitly re- 
solved [^,0. The pair potential is modeled by an ul- 
trasoft interaction which is logarithmic for an inner core 
and shows a Yukawa-type exponential decay at larger 
distances BUlT 



V ss (r) 
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-lnm + W- r<cr s ; 

W l + y/f/2 ~ S ' 

' cx P(~2^r( r _ else, 
(2) 

with k^T being the thermal energy. Our computer sim- 
ulations show that the so-called corona-diameter er s re- 
mains fixed for all considered arm numbers /, being re- 
lated to the diameter of gyration through er s ~ 0.66cr g 



f. However, the theoretical approach giving rise to eq. 
does not hold for arm numbers / < 10, because the 
Daoud-Cotton model of a star p3], on which the Yukawa 
decay rests, is not valid for small /. In these cases, the 
interaction inclines to a shorter-ranged decay for r > <r s . 
This is consistent with approaches in which at the linear 
polymer limit a Gaussian behavior of the pair potential 
is predicted (?]||^l| . Only the large distance decay of the 
interaction is affected; its form at close approaches has to 
remain logarithmic |22] ] . Accordingly, we propose the fol- 
lowing star-star pair potential for arm numbers / < 10, 
replacing the Yukawa by a Gaussian decay: 



V ss (r) = ~^k B Tp 



27^ ex P (-T 2 {r 2 - <7 S 2 )) else, 



(3) 



where t(/) is a free parameter of the order of l/i? g and 
is obtained by fitting to computer simulation results, see 
Fig. and Table 0. Using ra s (f = 2) = 1.03 we ob- 
tain for the second virial coefficient of polymer solutions 
the value -B2/-Rg = 5.59, in agreement with the estimate 
5.5 < B 2 /Rl < 5.9 from RG and simulations ||. 
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FIG. 1. Effective force between two isolated star polymers 
for arm numbers / = 2, 6 at center-to-center distance r. The 
theoretical result (line) derived from eq. (^) is compared to 
computer simulation data (symbols) to obtain the decay pa- 
rameter t. Rd is the nonvanishing core radius of one simu- 
lated star, values in Table [n|. Error bars are shown for the 
case / = 6 and provide an estimate for all /-values. 

An analytic form for the star polymer-colloid pair po- 
tential can be found by integrating the osmotic pressure 
of one star along the spherical surface of a colloid, fol- 
lowing an idea put forward by Pincus |23[ |. This can be 
achieved for arbitrary curvatures of the colloid but the 
analytical result below is accurate for size ratios q < 0.7 
and reads Hi : 
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where z — r — a c /2 is the distance from the center of the 
star polymer to the surface of the colloid. The constants 
are = 1/(1 + 2k 2 **) and £ 2 = & exp( K V s 2 )(l - 
erf(K(T s )). A(/) and «;(/) are fit parameters, obtained 
from computer simulations where the force between an 
isolated star and a hard flat wall is calculated, see Fig. ||. 
k is in order of l/<r g , see the values in Table |J, whereas 
geometrical arguments yield a limit Aoc = 5/36 for very 
large /. 



where g(x c , P, T) is the Gibbs free energy G(x c , N, P, T) 
per particle and P denotes the pressure of the mixture. 
If the function g(x c ) has concave parts the system shows 
phase coexistence and the phase boundaries can be calcu- 
lated by the common tangent construction on the g(x c ) 
vs. x c curves at constant P and T, ensuring the equality 
of the chemical potentials for each component in the two 
phases. The results obtained are shown in Figs. and H. 
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FIG. 2. Effective force between an isolated star polymer 
and a hard flat wall (q = 0) for arm numbers f = 2,6, 16 and 
32. z is the distance from the star center to the surface of 
the wall. Theoretical curves from eq. (Q) were compared to 
computer simulation data (symbols) to obtain the prefactor A 
and the decay parameter k. For better comparison we divided 
the force by f 3/2 . 

To access the thermodynamics of the mixture, we solve 
the two-component Rogers- Young (RY) closure, which 
is reliable for the one component star polymer system 
pB| and shows spinodal instability in highly asymmet- 
ric hard sphere mixtures |2(J. We also performed stan- 
dard Monte Carlo (MC) simulations using the interac- 
tions ([l])-(f§) as inputs and measuring the structure fac- 
tors at selected thermodynamics points. Excellent agree- 
ment between RY and MC was found. The thermody- 
namic consistency of the RY closure is enforced with a 
single adjustable parameter £; a simple scaling of the 
form £ Q/ 3 = £/o~ a p, (a, (3 = c, s) showed only small differ- 
ences compared to the unsealed form. 

The structure of binary mixtures is described by three 
partial static structure factors S a p(k), with a, (3 = c,s 
at wavevector k, obtained from the integral equations. 
Indication of a demixing transition is the divergence of 
all structure factors at the long wavelength limit k — > 
0, marking the spinodal line of the system. It is more 
convenient to consider the concentration structure factor 
Scon(k) = x 2 S cc (k) + x 2 S ss (k) - 2x c x s S cs (k). We denote 
the partial concentrations by Xi = Ni/N, (i = c,s). The 
k — > limit provides the approach to thermodynamics, 
through the relation : 
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FIG. 3. Binodals for the mixing-demixing transition in 
star polymer-colloid mixtures for different arm numbers 
/ = 2, 6, 16, 32 (from top to bottom) and size ratio q w 0.49. 
Symbols mark experimental results compared with theory 
(lines) for q = 0.50. 

Inside the spinodal line, the limits S a ,p{k —> 0) attain 
unphysical, negative values associated with the physi- 
cal instability of the mixture against phase separation. 
Consequently, a solution of the integral equations is 
not possible there, and above the critical pressure P* , 
S con (x c ,k = 0) is unknown in some interval Ax c (P). 
Thus, it is necessary to interpolate S con (x c ,k — 0) in 
order to perform the integration of eq. (|5|). In the vicin- 
ity of the critical point 77* ~ 0.3 the missing interval Ax c 
is very small and the interpolation is reliable. Here the 
binodals should be accurate, while for higher pressures 
(packing fractions rj c < n* and rj c > n*) the binodals are 
more approximate but show reasonable behavior. For 
highly asymmetric systems (q ^ 0.18) it becomes more 
and more difficult to get solutions of the integral equa- 
tions in the vicinity of the spinodal line and the calcula- 
tion of binodals was not possible. 

The results in Figs. || and [| show that theory and ex- 
periment are in good agreement. In particular, the same 
trends are found as functions of the system parameters 
/ and q. By increasing / and keeping q fixed, the demix- 
ing transition moves to lower star packing fractions rj s , 
as shown in Fig. 0. The largest differences occur at low 
arm numbers / < 10 which is theoretcially caused by the 
major changes in the star-colloid pair potential for low /. 
When q is decreased but / remains fixed, again a motion 
of the binodals to lower 77s is observed, as seen in Fig. 
|]. This trend is opposite to the one predicted by the AO 
model pS|]. The crucial parameter that determines the 
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trends of the phase diagram is the non-additivity of the 
mixture. Systems interpolating between the fully addi- 
tive hard sphere mixture and the fully nonadditive AO 
model can feature a depletion interaction which is more 
attractive than in the AO limit ^9|. The agreement be- 
tween theory and experiment is brought about without 
the use of any free parameters in the former that would 
allow for a rescaling of sizes or densities. All values are 
read off from experiment and the only free parameters 
of the theory appear on the level of the pair potentials 
and are used only in order to fit analytical expressions to 
the microscopically-determined star-star and star-colloid 
pair interactions. 
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T 




0.1 0.2 0.3 0.4 0.5 

Tic 

FIG. 4. Same as in Fig. ^| for an arm number / = 32 and 
different size ratios q. 

To summarize, we have derived from first principles 
the effective interactions between the components of a 
star polymer-colloid mixture, proposing analytical ex- 
pressions for these. We have used these to make theo- 
retical predictions about the thermodynamic behavior of 
the fluid phase of the system finding very good agreement 
with the experimental results we obtained independently. 
Further work regarding, e.g., the crystallization proper- 
ties of these mixtures is currently under way. 
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Monomer / M w ■ 10~ b [g/Mol] a) R c [nm] 



Rss 



la) 



PMMA 
PMMA 

PB 

PB 

PB 

PB 



104.0± 2.5 
289.0± 4.5 



16 

32 



0.57 ± 0.36 
1.51 ± 0.06 
3.45 ± 0.27 
4.87 ± 0.39 



51.0 ± 3.5 

52.1 ± 0.6 
51.1 ± 0.5 
51.4 ± 0.5 



a ) small angle neutron scattering (SANS) 
b ' static light scattering (SLS) 

TABLE I. Molecular characteristics of PMMA particles 
and PB star polymers. 
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A(/) 




ra s 


Rd/o- g 


2 


0.46 


0.58 


1.03 


0.04 


6 


0.35 


0.71 


1.16 


0.03 


16 


0.28 


0.76 




0.04 


32 


0.24 


0.84 




0.06 



TABLE II. Fit parameters A, k, t for the effective star-wall 
interaction of eq. (^) and the star-star interaction eq. (^|) ob- 
tained from molecular simulation. Rd is the nonvanishing core 
radius of one simulated star. 
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